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Abstract 

A full calculation of lepton-pair angular characteristics is carried out for e + e~ 
pairs created in p + Be and Ca + Ca collisions from 1.0 to 2.1 GeV/A. It is 
demonstrated that the dilepton decay anisotropy depends sensitively on the dif- 
ferent sources and may be used for their disentangling. Due to the dominance of 
the 77-and A-Dalitz decays and only a small anisotropy coefficient for it + it~ annihi- 
lation, the expected anisotropy coefficients show a decrease with invariant mass of 
the dilepton pair and change only moderately when comparing p + Be and Ca + Ca 
reactions at the same bombarding energy. 
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Dileptons are quite attractive electromagnetic signals since they provide almost direct 
information on the hot and dense nuclear phase in heavy- ion collisions at BEVAL AC/SIS 
and SPS energies Jj], 0, [|, f|]. The information carried out by leptons may tell us not 
only about the interaction dynamics of colliding nuclei, but also on properties of hadrons 
in the nuclear environment or on a possible phase transition of hadrons into a quark- 
gluon plasma (cf. ||). However, there are a lot of hadronic sources for dileptons because 
the electromagnetic field couples to all charges and magnetic moments. In particular, in 
hadron-hadron collisions, the e + e~ pairs are created due to the electromagnetic decay of 
time-like virtual photons which can result from the bremsstrahlung process or from the 
decay of baryonic and mesonic resonances including the direct conversion of vector mesons 
into virtual photons in accordance with the vector dominance hypothesis. In the nuclear 
medium, the properties of these sources may be modified and it is thus very desirable to 
have experimental observables which allow to disentangle the various channels of dilepton 
production. 

Recently, we have proposed to use lepton pair angular distributions for a distinction 
between different sources || [7| . Indeed, it could be shown that due to the spin alignment 
of the virtual photon and the spins of the colliding or decaying hadrons, this lepton decay 
anisotropy turns out to be quite sensitive to the specific production channel (cf. also HH). 
Whereas in our previous work || |7|, we have considered the dilepton anisotropy for 
elementary nucleon-nucleon or pd collisions, we now evaluate this new observable for the 
first time for proton-nucleus and nucleus-nucleus collisions from 1 to 2 GeV/A bombarding 
energy, where the calculated inclusive dilepton spectra can be controlled in comparison 
to the DLS data |,|,||. 

The dynamical evolution of the proton-nucleus or nucleus-nucleus collision is described 
by a transport equation of the Boltzmann-Uehling-Uhlenbeck type evolving phase-space 
distribution functions for nucleons, A's, iV*(1440)'s, iV*(1535)'s, pions and 7?'s with their 



isospin degrees of freedom. The details of the model are discussed in Refs. jry, [12]; here 
we use the same prescriptions and include additionally the production channels from the 
direct decay of the vector mesons p, u, and $ as well as the Dalitz-decay of the tu-meson 
|T3| employing the formfactors from Landsberg |I4"|. 

In Fig. m we present the calculated dilepton invariant mass spectra da/dM for p + 9 Be 
and 40 Ca + m Ca collisions at bombarding energies from 1 to 2.1 GeV/A and compare 
them with the experimental data of the DLS collaboration Jl|, [2|, |3| including the DLS 
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acceptance filter as well a mass resolution of 50 MeV. In the 'cocktail' plot the A Dalitz 
decay (labeled by 'A'), N* Dalitz decay ('JV*'), proton-neutron bremsstrahlung ('pn'), 
ttN bremsstrahlung (VAT') (all without VDM formfactor), 10 — > ne + e~ Dalitz decay, rj 
Dalitz decay ('77'), the pion annihilation channel (V7r') as well as the direct decay of the 
vector mesons ('p, u, $') are separated explicitely while the solid curves (denoted by ) aW) 
represent the sum of all sources. As seen from Fig. |l] the dilepton cross sections agree 
reasonably well with the experimental data as in |12| . We will use the same cross sections, 
however, without experimental DLS filter for the calculation of the anisotropy coefficients. 

In order to characterize the dilepton decay anisotropy we have introduced in || [7| the 
anisotropy coefficient B which allows to characterize the angular distribution of dileptons 
created in a hadron-hadron (h + h) reaction via 

da hh 

= rIMrl 1 ft = A ^ + ^ C ° S2 M- (1) 

dM a cos Vhh 

In Eq. (1) 9hh denotes the angle of the electron momentum l*_ - measured in the dilepton 
center-of-mass system (<f* = l*_ + — 0) - with respect to the velocity of the h + h 
c.m.s. (denoted by v hh ) relative to the dilepton c.m.s., i.e. v hh = —q h /qQ = —v^ h : 
cosdhh = {l-,Vg h ). Here, the momentum and energy of virtual photon (g^ 1 , (f 111 ) are 
defined in the c.m.s. of the h + h system while M is the invariant mass of a lepton pair, 
(M 2 = q 2 — (f). B\ h describes the anisotropy while A^ h determins magnitude of the cross 
section. 

The total differential cross section for h + h collisions now can be represented as a sum 
of the differential cross sections for all channels, 

Jf\ = E n = A hh (M)(l + B h \M) cos 2 M, (2) 

dMd cos dhh • , dMd cos hh 

rut i=criannel ™" 

which leads to the total anisotropy coefficient: 

da^ h B? h 



dM Til B hh 

B hh (M)= <B} h (M)>> <Bt h {M)>= , hh V , (3) 

i=channel flCr i _ 

^ dM ' 1 + l B hh 

where the special weighting factors originate from the necessary angle-integrations. Thus, 
the anisotropy coefficient B AB for A + B reactions is the sum of the "weighted" anisotropy 
coefficients (< Bf B >) for each channel i obtained by means of the convolution of Bf B 
with the corresponding invariant mass distribution (cf. |]10[). 
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For heavy-ion reactions the situation becomes more complicated due to the nuclear 
dynamics and the explicit time evolution of the interacting system. Here we start from 
the point that the form of the angular distribution for all 'elementary' interactions a + b, 
that occur in the nucleus-nucleus reaction A + B, are known. For this aim we employ the 
results of our previous works || |7|, [L(J where the anisotropy coefficients for the pp, and pn 
bremsstrahlung, NN — > A — ► NNe + e~ Dalitz decay, rj Dalitz decay and annihilation 
channels were calculated explicitly in the hadron-hadron center-of-mass system a + h. 

The differential angular distribution in elementary a + b collisions - before averaging 
over the momentum g ttb - can be represented as: 

daf b 



A? (l+Bf cos 2 6 ab ), (4) 



dM dq ab d cos Q a 

where the coefficient Bf b for the elementary process a + b is a function of the dilepton 
mass M, the masses m a ,mb and the initial invariant energy y/s of the hadrons involved 
in the reaction: Bf b = B^ b (M, \fs]m a ,m^). It is important to note that we know the 
'elementary' coefficient only in the a + b system. In nucleus-nucleus collisions, however, 
the direction = q ab jq^ is changed in each elementary a + b collision. Thus in order 
to define an anisotropy coefficient in the latter situation we have to perform an angular 
transformation from the elementary c.m.s. (Q a b) to the c.m.s. of the colliding nuclei {9 ab)'- 

cos Qab = cos d ab cos 6 g + sin d ab sin 8 g COs((fAB - Vq) ) (5) 

where 9 q is the angle between the dilepton c.m.s. velocity v^ b = ^/q^ in the c.m.s. of 
a + b and the dilepton c.m.s. velocity v^ B = q AB /q^ B in the c.m.s. of the colliding nuclei 
A + B: cos^g = Vq b ■ v AB /(\Vq b \ \v£ B \). Here, the vector q^ b is obtained by a Lorentz 
transformation of q AB : q ab = L(V ab )q AB , where V ab = (p a + Pb) / (E a + E b ) is the velocity 
of the a + b system relative to the A + B system. 

Substituting (||) into Eq.(£|) and using dflAB = dVL ab we get the respective distribution 
in the observable angle Qab by integrating over the azimuthal angle <Pab'- 

daf 



dMd cos Qab 



- J dq AB A t (Q q ) (l + B t (Q q ) cos 2 Q AB ) , (6) 

Dab R a ^ 

MP*) = 1 + ^-sin 2 fl g , B(Q q ) = -^-(3 C os 2 Q q - 1). (7) 
Obviously, the transformation (^|) does not change the quadratic form in cos Q of the 
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angular distribution. From Eqs. ([!]) and (§) and the normalization condition we finally get 

^ = z^b ' ( 8 ) 



dn ab 

fdq AB A(6 q ) ' 



dMdq AB 

where the differential cross sections daf b / (dMdq AB ) for the different channels i are taken 
from the BUU calculations as described in [11], The respective 'reduced' differential 
cross sections do~ AB jdM are shown in Fig. [I] for the systems to be studied below. 

The above definition of the anisotropy coefficient in nucleus-nucleus collisions is valid 
for all channels except for 7r + 7r~-annihilation because there are only two particles in the 
initial and final states. For this particular reaction we have to use the angle of the pion 
momentum with respect to the lepton momentum in the c.m.s. of the leptons (or pions, 
which is the same for this channel), i.e. q ab = p a + p = l + + l_ = 0. As was shown in 
Ref. M, the elementary anisotropy coefficient in the tttt c.m.s. then is B n + n - = — 1. In 
heavy-ion collisions we can use the same definition for the polar angle as for the other 
channels because we can reconstruct the direction of q AB in the c.m.s. of the colliding 
nuclei A+B. Furthermore, we performed the same angular transformation from the c.m.s. 
of the pions a + b to the c.m.s. of the nuclei A + B as in Eq.([5]) with the replacement 
q — > n , Here, 6 n is the angle between the pion momentum p a in the c.m.s. of a + b and 
the vector v AB = q AB jq§ B ■ Following the same angular integration as for the previous 
cases we end up with the expression for the anisotropy coefficient for 7r + 7r~ annihilation: 

b ab = JdcosO^ A(6 W ) B(e w ) ly(cos^) 
JdcosO^ A(6 n ) W(co8 9 n ) 

The quantities A(9 n ), BiO^) are defined by Eq.(0) using B ^ = —1 while cos^ can be 
computed by the scalar products of the four momenta of pions p a ,Pb and colliding nuclei 

Pa, Pb, 

n (p a - Pb ).(P A + P B ) M 

COS0,r = i , (10) 



\ M 2 _ 4m 2) ( [(po +pb ).(p A + p B )f _ [P A + p B )2 M 2 

with M 2 = s a b = (p a + Pb) 2 following O.V. Teryaev [ 15[ . 

A closer look at Eq.(|9p shows that B AB is a constant; its absolute value follows from 
the pion angular distribution W(cos8 n ). In case of an isotropic angular distribution, 
^(cos^) = const, the anisotropy coefficient B AB = according to Eq.(|9|). Moreover, in 
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the BUU calculation the angular distribution W(cos6 n ) is also a function of time t due 
to the dynamical evolution of the system. Technically we define the distribution W(t) 
at time t (i.e. in the time interval [t — At/2 : t + At/2]) as the ratio of the number of 
pions with fixed cos 6 n produced in the above time interval to the total number of pions 
produced in the A + B reaction, i.e. W(t, cos9 n ) = N(t, cos9 v ) /N tot . 

In order to demonstrate the pion angular anisotropy we show in Fig. |^ (l.h.s.) the 
pion angular distribution W(t, cos 9 V ) for p + 9 Be at 2.1 GeV and 40 Ca + m Ca at 1.0 and 
2.0 GeV/A at those times t when the maximal number of pions N(t) was produced during 
the interval At; the respective 7r-production rate N(t)/N to t is illustrated in the r.h.s. of 
Fig. |^ for the same reactions. Due to the low number of pions produced for p + 9 Be 
the angular distribution suffers from low statistics and the error bars - which result from 
different runs with 8000 testparticles per nucleon - are large. For the further analysis 
we have used the dotted line in Fig. |2| which represents a fit to the angular distribution. 
On the other hand, for 40 Ca + 40 Ca collisions the statistics is good enough and we show 
only fits to the 'numerical' data in terms of the solid and dash-dotted lines, respectively. 
The anisotropy of W(cos9 n ) is most pronounced for Ca + Ca at 1 GeV/A, because the 
probability for 7r's with transverse momentum is larger then with longitudinal momentum 
due to a stronger pion absorption in beam (longitudinal) direction. At 2 GeV/A, the pion 
distribution becomes more isotropic in line with a 'pionic fireball' scenario. Since pion 
absorption effects are only small for p + Be, the resulting angular distribution shows only 
a very modest anisotropy. 

Before going over to the calculation of the anisotropy coefficient for nucleus-nucleus 
collisions one has to take into account that resonances (A, N*) can be created in quite 
different elementary channels than in pN or pd reactions. For example, A production 
via the 7rN — > A channel becomes quite important; an elementary channel for which the 
dilepton anisotropy has not been computed so far. We thus have calculated the e + e~ 
anisotropy using the same vertices, delta-propagator, coupling constants and formfactors 



as in Refs. 0, [Lb]. Since this evaluation is straight forward, we do not present the details 
here. 

The results of our microscopic calculation for the anisotropy coefficient for the 7r+ N — > 
A — > e + e~ + X channel are displayed in Fig. ||. The anisotropy B nN _>& is a function 
of the dilepton invariant mass M and the invariant energy of the interacting particles 
s = (p w +pn) 2 = M\. For small yfs only deltas with M A rs M ao = 1.232 GeV appear 
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and the coefficient B n N->A ^ 1 as expected for the Dalitz decay of a free delta || . With 
increasing more energetic deltas can be created and (for fixed M) the phase-space for 
the final nucleon and virtual photon increases leading to a decrease of B^n^a- We note 
that we do not take into account the Dalitz decay of the higher nucleon resonances in the 
further calculations because their statistical weight is too low. 

Fig. f| finally shows the computed weighted anisotropy coefficients < Bi(M) > for 
p + 9 Be and m Ca + i9 Ca collisions at the bombarding energies from 1 to 2.1 GeV/A. The 
main contributions arise from the 77 and A Dalitz decays due to their large 'elementary' 
anisotropy coefficients and cross sections (cf. Fig. [I]), respectively. The contribution 
from pn bremsstrahlung is practically zero at all energies due to a smaller 'elementary' 
anisotropy coefficient and due to a lower cross section as well. The weighted coefficient 
from 7r + 7r~ annihilation is rather small (~ 0.1) even for the Ca+Ca reactions and decreases 
for M > m p due to the threshold behaviour of the cross section. However, compared to 
< B n + n - (M) > for p + 9 Be, where pion annihilation is very low, a clear (but moderate) 
enhancement can be extracted. 

The contributions of the further channels (N*, ttN bremsstrahlung, u — > 7r°e + e~) are 
also negligible due to their small cross sections (cf. Fig. [I]). The cross section from direct 
decays of the vector mesons becomes compatible with pion annihilation for M « m p for 
p+ 9 Be at 2.1 GeV, but the anisotropy coefficient for the 'free' vector meson in the vacuum 
is zero [17]. We do not discuss here a possible modification of the p-meson properties in 
the medium that might lead to non-isotropic angular distributions of dileptons because 
for p + 9 Be at 2.1 GeV one cannot reach sufficiently high baryon or pion densities. 

On the other hand, for 40 Ca + 40 Ca collisions at 2.0 GeV and for p + 9 Be at 2.1 GeV 
we observe a clear anisotropy from the rj Dalitz decay, while for p + 9 Be at 1 GeV the A 
Dalitz decay gives the main contribution at small invariant mass due to a dominant A 
cross section; the 77 coefficient increases at M from 0.4 to 0.5 GeV for the same reason. 

For clarity we briefly discuss the steps to 'extract' the anisotropy coefficient from 
experimental data. In dilepton experiments the four-momenta of leptons in the c.m.s. 
of the nuclei A + B (or in the laboratory frame, which is connected with the c.m.s. 
by a simple Lorentz transformation) are measured - (e AB ,l AB ), (e AB ,l AB ). The four 
momentum q = (q AB , q AB ) of the virtual photon in the A + B c.m.s. then is given by: 
q^B — £ ab _|_ £ ab^ qAB _ j^ab _|_ rj^g arL gi e q ab t Q be computed via: 

\q AB \ -q£ B cos9 ql 
C ° SdAB = \qn C os9 gl -q^ (H) 
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where 9 q i is the angle between the q AB and l AB , i.e. cos9 q i = l AB ■ q / \ \q AB \). 

In order to extract the anisotropy coefficient one has to count the number of dilepton 
events N(M, cos 6^) with fixed cos8ab and invariant mass M, The anisotropy coefficient 
for the invariant mass M then is simply given by: 

B {M) - N(M,cos6 AB = 0)- h (12) 

Thus, summarizing, the calculated anisotropy coefficients for p + Be and Ca + Ca 
collisions support our suggestion in Refs. [|[ [7|, [H| that the dilepton decay anisotropy 
may serve as an additional observable to decompose the dilepton spectra into the various 
sources. Since the anisotropy vanishes in a hadronic fireball scenario our present results 
provide valuable information about the nonequilibrium stage of the reactions. 

We gratefully acknowledge many helpful discussions with Gy. Wolf, S.S. Shimanskij, 
S. Teis, O.V. Teryaev, A.I. Titov, and D.V. Toneev. 
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Figure captions 



Figure 1: The dilepton invariant mass spectra da/dM from p + 9 Be and m Ca + 40 C a 
collisions at bombarding energies from 1 to 2.1 GeV/A in comparison to the experimen- 
tal data [III |2|, The "77" denotes the contribution of the //-channel, the "A" labels 
the contribution of the A Dalitz decay, "iV*" the N* Dalitz decay, u uj — > 7re + e~" is 
the uj — > 7re + e~ Dalitz decay, "pn" the proton-neutron bremsstrahlung, u irN" the 7riV 
bremsstrahlung, "7T7r" the pion annihilation channel, while u p,u>, $" is the direct decay of 
vector mesons . The solid curves (denoted by "a//" ) show the sum of all sources. 

Figure 2: The pion angular distribution W(t, cos6 n ) for p + 9 Be at 2.1 GeV and m Ca + 
40 Ca at 1.0, 2.0 GeV and time t (l.h.s.) where the maximal number of pions N(t) was 
produced. The pion production rate N(t)/N to t is shown in the r.h.s. of the figure. 

The results of our calculation for the anisotropy coefficient of the tc + iV — > 
- + X channel at M A = y/s from 1.232 to 1.8 GeV. 

The weighted anisotropy coefficients < B^M) > for p + 9 Be and 40 Ca + 40 C a 
at bombarding energies from 1 to 2.1 GeV/A. The notation is the same as in 



10 



0.020 



0.015 - 



<x> 

© 0.010 



t=7 fm/c, p+ v Be, 2.1 GeV 



__t=9.5 fm/c, 40 Ca+ 40 Ca, 2.0 GeV 



0.005 - 




0.000 



-1.0 



■0.5 



0.0 

cos 9 



0.5 



1.0 



71 



Fig. 2 





M (GeV) 

Fig. 3 



